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ABSTRACT 


The fracture problem of laminated plates which consist 
of bonded orthotropic layers is studied. It is assumed 
that the medium contains periodic cracks normal to the bi- 
material interfaces and the external loads are applied away 
from the crack region. The fields equations for an elastic 
orthotropic body are transformed to give the displacement 
and stress expressions for each layer or strip. The un- 
known functions in these expressions are found by satisfy- 
ing tne remaining .boundary and continuity conditions. A 
system of singular integral equations is obtained from the 
mixed boundary conditions. Three cases are considered: 

a) The case of internal cracks 

b) The case of broken laminates 

c) The case of a crack crossing the interface. 

The singular behavior around the crack tip and at the bi- 
material interface is studied. It is shown that the crack 
surface displacement derivative has a power singularity for 
practical orthotropic materials when the crack touches the 
interface, i.e., for case (b). In studying the singular 
behavior at the bimaterial interfaces in case (c), it is 
found that for some orthotropic material combinations there 
is no singularity in the crack surface displacement deriva- 
tives and the stresses. In each case the stress intensity 
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factors are computed for various material combinations 
and various crack geometries. The results for orthotropic 
materials are discussed and are compared with those for 
isotropic materials, 



1. INTRODUCTION 


In structural design one of the most important consid- 
erations is the fracture of individual components. Al- 
though, fracture may not always mean total failure, it is 
considered in modern engineering as an important problem 
for safe and economic design of structures. It would be 
very attractive to develop special types of designs for 
which the structural resistance to fatigue crack propaga- 
tion is improved. In the aerospace industry, the use of 
composite sheet materials with buffer strips parallel to 
the main load-carrying laminates seems to be such a design 
practice. The process of manufacturing composites gives 
the opportunity to improve the structural resistance to 
fatigue crack propagation by strengthening the material in 
certain directions. The increasing use of composites in 
structures generates new problems for the structural de- 
signer. Among these problems, we are mainly interested in 
the fracture of layered composite materials. 

There are two main problems in studying the fractpre 
of composites: the development of an appropriate failure 

criterion and a mathematical model for the calculation of 
the related load factor. The failure criterion affects 
the course of the analytical work in the sense that it is 
the failure criterion which generally determines the phys- 
ical quantities that one should compute (such as the stress 
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intensity factor, the strain energy release rate, COD, 
etc.)* There are many failure criteria or theories which 
are used to predict failure of structures. In Elastic 
Fracture Mechanics where only small scale yielding is al- 
lowed, K < Kjq is such a criterion. In this case failure 
occurs when the calculated value of the stress intensity 
factor reaches a critical value, Kjq, which can be deter- 
mined experimentally as a material property. There are 
also other one-parameter failure criteria (such as critical 
plastic stress intensity factor Kp^ and J integral) which 
have been recently proposed to predict failure from elas- 
tic to fully plastic range. K is a very highly effective 
correlation parameter in studying the fatigue crack prop- 
agation phenomena. In aerospace structures the basic 
problem is the nucleation and propagation of fatigue crack 
which may eventually reach a critical size causing cata- 
strophic failure. That is why, in this study we focus our 
interest to the computation of the stress intensity factors 
and in the investigation of the singular behavior of the 
stress state around the crack tips. 

In studying the fracture problem of composites, a 
mathematical model, which will reflect the geometrical and 
physical properties of the medium and the real mechanism 
of fracture, is needed. Because of mathematical diffi- 
culties and the lengthy computation that the analysis 
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requires, in the recent studies the geometry and the ma- 
terial properties have been considerably simplified. The 
problem of a multi-layered isotropic medium, which corMsts 
of many layers and where a crack normal to the interface 
can appear, has been treated by Hilton and Sih [1]. In 
this problem the geometry is simplified to a single layer 
between two dissimilar half-planes where the elastic prop- 
erties are averaged. The same problem has been considered 
by Bogy [2]. The problem of a broken laminate between two 
half planes has been investigated by Ashbaugh [3] and Gupta 
[4]. The extension of the problem treated by Hilton and 
Sih to orthotropic media has been solved by Arin [5]. The 
fracture problem of a composite plate which consists of 
parallel load-carrying laminates and buffer strips has re- 
cently been solved by Erdogan and Bakioglu [6]. In this 
work the load carrying laminates and buffers are consid- 
ered to be isotropic and linearly elastic. The orthotropic 
case of the problem treated in [3] and [4] has also been 
solved by. Arin [7]. 

The objective of this work is to investigate the 
fracture problem of composite plates containing periodic 
buffer strips. The laminates and buffer strips are assumed 
to be linearly elastic and orthotropic. In general, this 
is the case in the actual plate and shell structures such 
as those, for example, which consist of boron-epoxy 
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composites. It is also assumed that the fatigue cracks 
may appear and propagate in main laminates, in buffer 
strips or in both normal to the interfaces. The external 
load is applied to the plate parallel to the strips and 
away from the crack region. Three different problems are 
studied: the internal crack problem, the case of broken 

laminates and the case of a crack crossing the interface. 
A general formulation of the problem is given for plane 
strain and generalized plane stress cases by the use of 
Fourier Integral Transform Technique. The singular be- 
havior around ends and at the bimaterial interfaces is 
studied. The resulting singular integral equations are 
solved numerically and the stress intensity factors are 
calculated for various crack geometries and various ma- 
terial combinations. 








2.4 The Field 


For an orthotropic solid the matrix [A-jj] is: 


r ” 


11 

Ai2 

^13 

0 

0 

12 

A22 

^23 

0 

0 

1 3 

^23 

^33 

0 

0 

0 

0 

0 

A44 

0 

0 

0 

0 

0 

^55 

0 

0 

0 

0 

0 


Defining the inverse of [A^-j] by 


[aij] = [Aij]' 


for orthotropic materials we have: 

1 _ ^zx 

^ ^y" ^z 

. ^xy _L ^zy 

£ J E w E 7 


V V 

_ xz ^ yz 


0 0 0 

0 0 0 

0 0 0 

1 rt n 


XZ 

0 0 


^x^yx ' ^y^xy * ^v^zy " ^z^yz ’ ^z^xz ^x^; 


y zy 


( 2 . 5 ) 
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Substituting (2.2) into (2.3) and using (2.1) and (2.4), 
the stresses and the equilibrium equations can be expressed 


in terms of the displacements as follows; 


A M+a ^ 

''ll 3x ''iz 3y ''is 3z 


A £ii+A ~+A — 

''12 3x ''22 3y ^ ''23 3z 


3u 

3x 

3U 


A ^+A — +A — 

''13 3x ''23 3y ^33 3z 


■yz 


A (^1 + 
''44^3z 3y' 


xz 


A + M) 

''55'3x 3z' 


xy 


A (M + IX) 
''66^3y 3x^ 


(2.6) 


3^u 


3^u 


3^u 


+ (A-|«+Acc) + (A-l;,+Arc) 


3^w 


^66 w ^ ^5 Ip- ^ '”12""66^ W " '"13^"55^ ^ 


A liv . llv + A + YA +A 1 + f A +A ) 

''66 ^ ''22 3y^ ''44 ^ '''l2 ^'66^ 3x3y ^''23 ''44^ 3y3z 


0 

0 


I 3^W 

55 W 


+ A 


3" 

44 3? 


w 


+ A 


3^w 
33 ^ 


^^3‘^^55^ 


3^1 


3x3z 


(A 23 +A 44 ) 


3y3z 


0 


(2.7) 


2.4.1 Case of Plane Strain 

For the plane strain case we have: 

u = u(x,y) , V = v(x,y) , w = 0 (2.8) 

and from (2.2), 
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(2.9) 


9U _ 3v _ 9u . 9V 

" 3x ’ ~ 3y * ^xy ~ 3y 3x 


^ ’ ^yz ^ ’ *^xz ~ ^ 


Thus, tha stress-displacement relations and the equilib 
rium equations become: 


^x 

< 

II 

M + A 

9x ^12 

3v 

3y 


"y 

= A^2 

3x ^22 

3v 

3y 


^z 

" ^13 

^ + A 
3x ^23 

3v 

ay 



^66 

/ 3u . 3Vx 
' 3y 3x^ 



^yz 

^xz 

= 0 


(2 



3^U . 
3x^ 

*66 

3^U , 

* 

(Ai2 


3^v 

3x3y 

*66 


*22 

+ 

^^1 2 


3^u 


3y^ 

3x3y 


2.4.2 Case of General 1 zed Plane Stress 

In this case since a = x _ = x = 0 from (2.3), 

Z /\ Z jr iu 

the average stresses and strains we can write: 


r 1 


- - 

e 



X 


x 


= [a] 

'’y 

Y 

’xy 


''xy 

_ _ 


— 


. 10 ) 


. 11 ) 


for 
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where 


[ 5 ] = 


1 

r 


X 

0 


0 

0 

1 


. [A] = 


1 

V 

M. 


0 


E «A 
X 


1 

0 


0 


xy 


and A = 


1 


“ ^xy^yx^- 

X y 


(2.12) 


The equilibrium equations reduce to: 

N 

3^v 
33^ 3x3y 

3^u 


^1 

3^u . 
3)F ^ 

^33 

3^U . 

w 

(Ai2 

■*■ ^33 

*33 

3^v . 
Tx^ 

A22 

3^V . 

Jy^ 

(Ai2 

^33 


= 0 


(2.13) 


Considering the structure of equations (2.11) and (2.13), 
the equilibrium equations can be written for plane strain 
and generalized plane stress cases in the following form: 

, i^u ^ l^u ^ fl ,9X = 


OF 


3^U . 3^U . o g V _ A 

+ ly^ ^ ^3 3y3X ^ 


a^'v . A Xv + R 9^^ = 0 

■^ + 32 ayT + 33 axay 


(2.14) 


where 


3 

and 

3 


1 


ii 

66 


'll 


1 


A 33 


, 3 o “ 


^22 ^ _ 1 + ^^12. for plane strain 


. 3> 


^•22 

A33 


At o 


» 3 a = s 


A 33 


for plane stress. 

(2.15) 
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3. 


DISPLACEMENT AND STRESS FIELDS FOR STRIPS 


The two-dimensional composite medium is formed of two 
sets of periodically arranged strips having widths 2h^ and 
2hg as shown in Figure 1. They are perfectly bonded along 
their straight boundaries, and contain symmetrically lo- 
cated cracks normal to the interfaces, of length 2a and 2b 
respectively. The load is applied away from the crack 
region, such that the crack plane is a plane of symmetry. 

Using the usual superposition technique, the solution 
of the actual traction-free crack problem may be obtained 
by superposing the homogeneous uncracked strip solution to 
the solution of a cracked strip loaded with self-equilibra- 
ting crack surface tractions (see Figure 2). Since we are 
interested only in the computation of stress intensity 
factors and the singular behavior of the stresses around 
crack ends, we will consider only the singular part of the 
solution, where the self-equilibrating crack tractions are 
the only external forces. 

First we will find solutions to (2.14) satisfying 
certain boundary conditions of the strips. The combination 
of these sol utions wi 1 1 be forced to satisfy the remaining 
boundary and continuity conditions. 
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3.1 Solutions u^^^(x,y), v^^^(x,y) 


Assume: 

u^®hx,y) = I f f(a,x)cos ay 
J 0 

v^^Nx.y) = - f g{a,x)sin ay 
1 0 

Substituting (3.1) into (2.14) we obtain: 

df _ 


da 

da 






The solution of (3.2) can be written as: 


f(a,x) = A(a)e^'“^ + B(a)e’^^"^ + C(a)e^ 20 tx ^ D(a)e 


s lax 


g(a,x) = 3y A(a)e^^“^ - B(a)e 


- D(a)e'"^“^ 


where Si and S 2 are the roots of 


s- + B4S2 + 35 = 0 


+ 6 , 


C(a)e 


S 2CtX 


and 


3] 32 


/3^ - 43'i 

4 


^3^1 


1 - 3iS 


1"2 


33S2 
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(3.1) 


(3.2) 


Siax 


(3.3) 


(3.4) 



From (3.4) we can write 


= w^ + iWg = V(- 34 + 3 g )72 


So = + iw. = /(- 3/, - 


S3 = > S) . 


s-j and $2 are both real or complex conjugates 


3.2 Solutions u^^^(x.y). v^*^Mx.y) 


(b) 


(3.5) 


Assume: 


u<‘’>(x,y) 


v<‘’>(x.y) 


2 
7T 

2 ( 

IT 


/■oo 


h{a,y)sin ax da 


S,{a,y)cos ax da 


Substituting (3.6) into (2.14) we have: 


(3.6) 


d^h 

dy2 


3-5 a 


dil 


^ " ‘=1 


3,a^ h = 0 


o d^iZ> .0 ^ dh ^2n - n 

®2 3y^'^B3“3y'“‘‘ 

Solving (3.7) we obtain: 

h(a,y) = E(a)e^'“y/'^ + F(a)e‘^ ^ 

+ 6(a)e"^^y/’^ + H(a)e-^^'^y/*^ , 


(3.7) 


and 


5^(a,y) = 3g 


•^(^)^s,ay//67 . r(a)e-"-“y''*^l 
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where 




^9 ° b\ 
- A. 

10 “ 3a 


3i/3s ' 

^2 /F7 


(3.8) 


A superscript * will be used for the material constants and 
unknown functions when the above expressions are used for 
the second strip. 


If one examines the roots of equation (3.4), he will 
realize that there are two types of orthotropic materials. 
We will denote the material as type I when s^ and S 2 are 
real, and as type II when they are complex conjugates. We 
will assume, in our analysis, that the materials of both 
strips are of type I. Similar analysis can be done for 
the remaining combinations. 

3.3 The Displacements 

For each strip, we can write: 

u(x,y) = u^^^(x,y) + u^^^(x,y) 

v(x,y) = v^®^(x,y) + v^^^(x,y) . 

Noting that; 

u(x,y) = - u(-x,y) and v(x,y) = - v(x,-y) 
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We will obtain: 


B{a) = - A (a), D(ci) = -C(a), F(a) = -E{a), H(a) = - G(a). 


For material type I ; 


^1 ~ ~ 3nd real. 


g,P7,M3.Mg,HlO 


Using the information given above, and keeping in mind that 
u and V vanish when y goes to infinity, for y>0, the dis- 
placement expressions can be written as follows: 

-oo 

u(x,y) = ^ j [A(a)sinh(w^ax) + C(a)sinh(WjOx)]cos ay da 


+ — 
IT 


0 


sin ax da 


v(x.y) - i 


^[3yA(a)cosh(Wiax) + 3 gC(a)cosh(w 2 ax)]sin ay da 

sign(wi )3gE(a)e“^'^^ + sign(w2)G(a)g^^l'^^^°^'^^*^ 

•cos ax da (3.9) 


3 . 4 The stresses 

For generalized plane stress case: 

^11 ^x ^ 2^ y 

^xy " ^33'^xy ' (3.10) 
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Differentiating (3.9) and using (2.12) and (3.10), 


the stress expressions can be written as: 




y^E{a)i 


Iwilay//^ ^ 




Iwalay/i/^ acos ax da + |[ 2 Y 3 A(a)cosh{Wiax) + 2v^C(a) 


•cosh(w 2 ax)]acos ay da 


tt(1-v V ) 




ay(x.y) = I Y5E(a)e-l'^^l«y/»^ + YgG(a) . 


CO 

acos ax da + j [2YyA(a)cosh(w-|ax) + 2YgC(a) 


cosh(wjax)]acos ay da 


Tj^y(x,y) = [2Y9A{a)sinh(w-,ax) + 2 y^o^^“^ 


sinh(w 2 ax)]asin ay da + 


'11 






asin ax da 


(3.11) 


These expressions are valid also for the plane strain case 
with the following substitutions: 

''yx ' ’ ''xy ° = I/A-,, . 

(Ej,-A) = VA22 . 

The elastic material constants yj ^^e defined in Appendix A. 
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FORMULATION OF THE PROBLEM 

The solution of the problem may be obtained by de- 
termining the unknown functions which appear in the dis- 
placement and stress expressions, under the following 
boundary and continuity conditions: 

u^ (h^ ,y) = u 2 (-h 2 ,y) , 

v^(h^,y) = V2(-h2»y)(0 < y < “) (4.1a,b) 

'^Ixy^^l*^^ ^ '^2xy^"^2*^^ ^ ^ (4.2a,b) 

u^(0,y) = 0 , T^^y(O.y) = 0 (0 < y < co) ( 4 . 3a , b ) 

® » ’’■2xy^*^’^^ = 0 (0 < y < «) (4.4a ,b) 
'T-j xy ^ ® ~ 0 , lxil<h.| , 

i 2 xy(x 2 . 0 ) = 0 . Ix 2 i<h 2 (4.5a,b) 

<^ly(Xi, 0 ) = -pi(x-j) , [x-jl<a , 

(x^ .0) = 0 , a<| x^ l<h^ (4.6a,b) 

V 2 (X 2 , 0 ) = 0 , b<lx 2 l<h 2 (4,7a.b) 

The conditions (4.3a,b) and (4.4a,b) are satisfied iden- 
tically. 
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Using (4.5a,b) we obtain: 

G(a) = - -il. E(a) , 6*(a) = - E*(a) 


1 2 


Yi2 


The mixed condition (4.6) gives: 


11m f E(a) 

A+ •'0 


Yce 


- wjay/i/e7.„ Lj- e"|w3|ay/»^ 


^6 


Yi2 


acos ax-jda 


I [2YyA(a)cosh(w-jax^ ) + 2 YgC{a)cosh(w 3 ax^ )]acos ay da 


• |Xll<P 


(4.8a) 


and 


o 

v-|(x^,0) = - ^Yi 3 E(a)cos ax^da - 0 a<lx^[<h^ 

(4.8b) 


Define , 


9Vi(xijO) 

^ = 4>(x-|) such that <f>(x^) = 0 for jx^|>a. 

(4.9) 


Differentiating (4,8b) with respect to x^ and taking the 
inverse transform, we obtain: 


Y-| 3«E(a) 



)s1n ax^ 



(4.10) 


If we now substitute (4.10) into (4.8a) and evaluate 
some of the integrals in closed form (see Appendix C) we 
will end up with the following singular integral equation: 
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14 ] dt + [2YyA(ot)cosh(w-jax^) + 2YgC(a)cosh(w2aK^ )]ada 


7T(1-V V ) 

- - — ,<a 

where because of symmetry (J>(t) = - (|)(-t). 


(4.11) 


Similarly defining. 


9V«(Xp,0) * * 

— 1^^ =<t>(x 2 ) such that(j>(x 2 ) = 0 for |x 2 l>b 

(4.12) 

and using the mixed condition (4.7) by the same procedure 
we obtain: 


dt + J [2YyA (a)cosh{w^ax2) + 2 yqC (a)cosh(w2aX2 ) ]otda 


^14 


-b 


tt{1-v* V* ) 

P3(x,) 


-b<X2<b 


(4.13) 


The next step is to determine the unknown functions A(a), 

‘At 'fc 

C(a), A (a), C (a). This can be done by using the contin- 
uity conditions (4.1a,b) and (4.2a,b) and taking the in- 
verse transforms. Then we obtain the following system of 
1 inear equations : 


2A(a) sinh(w^ah.| ) + 2C (a)sinh(w 2 ah^ ) 

'ff ■ ’fc 

+ 2A (a)sinh(w^ah2) + 2C (a)sinh(w3ah2) = Ri(a) 
2$yA(a) cosh(w-| ah-| ) + 23 gC(a)cosh(w 20 ihi ) 

ic if if if if if 

- 2^j A (a)cosh(w-|ah 2 ) - 23gC ( a) cosh ( W 2 ah 2 ) = R 2 (a) 
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2Y3A(a) cosh(w^ah^ ) + 2y^C(a) cosh( w^ah^ ) 




★ 'k 


- 2\^y2^ (a)cosh(w-jah2) - 2X^Y4C (a)cosh(w20th2) = 
2YgA(a)sinh(w^ah^ ) + 2Y^QC{a)sinh(w3ah^) 

4r 4r “k k k k 

+ Zy^X^P^ (a)sinh(w-jah2) + 2 y^qX 2C (a)sinh(w3ah2) = R^Ca) 

(4.14) 

R.(a) and X. are defined in Appendices B and A respectively. 

’ J 

Solving (4.14) we obtain: 

1 R 1 ( *^ ) Rp ( ^ ) 

" 2cosh{w^ah-, ) f(a) ^2''^^ ^ f(a) ^2^“^ 

Rsta) R4(‘^) 

f(a) "^2^“^ f(a) "2^“^ 


~ 2cosh(w2ah^ ) 
Ro(a) 


If > % 

A (a) = 


R-j { ct ) Rrt ( ct ) 

f{a)“ 9i(“) f(a) 


RgCa) 


2cosh(w*ah2) 

H 

C (ot) = ; 

2cosh(w3ah2) 

Rq(a) 


R4(a) 

■ 

(a) 

f(a) 

R,(a) 

9o(“) 


f (a) 


R4(“) 

TT5T 

•no(«) 

R,(a) 

g(a) 

+ 

TTO" 


R4(“) 

« 

n(a) 


f (ct) 



Ro{ot) 


Ro (a) 


(4.15) 
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The functions f(a), g(a) etc., used above are given in 
Appendix B. 

Substituting (4.15) into (4.11) and (4.13) we obtain 
the following system of singular integral equations: 



j k-|’j(xi,t)(J)(t)dt 


(t)dt = 2 y^Je^Pi^^1^ 


ib 


- a < x-| < a 


1 

IT 


rb */. \ /-a 

l-^dt + 
t-x, 

-b 


k2i(x2,t)4)(t)dt 

I k22(x2.t)({)*(t)dt 


-a 

b 


-b 


P2<*2) 

2v,4e; 


b < X 2 < b 


(4.16a,b) 


where 


kii(xi.t) = 


,00 


14 




^ k2(x 


da 


l<12(xpt) - ^ 


14 


k3(x, 




w„ 


da 
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fOO 




+ k5(X2.a)e-'">-t>“>^/l*'’l 


da 


k22(x2.t) = j“ k2(x2,a)e-(''2-‘)“'^/|“t 

14 ‘^L 

+ kg(x2,a)e-‘'’2-‘)“'^/l«? 


(4.17) 


The expressions k.(j =1,8) are given in Appendix B, 

J 

By letting h 2 ‘>'“ (or h-|->-“) one can recover the spec- 
ial case studied in [5]. For a < h-j and b < h 2 » the in- 
tegrands of kernels k..(i,j = 1,2) vanish when a -»■ “ and are 

• w 

bounded for all values of a, except when a = 0. 


Around a = 0 the asymptotic behavior of the inte- 
grands I., of the kernels k.. is of the following form: 

* si * J 

^ + 0(1) (l.j =1.2) (4.18) 

where the c.. 's are known constants. In order to obtain 

I J 

a solution, one should show that the singularity due to 
1/a is removable. Consider the following integral: 


ra 


r“ /-a 

(|>(t)dt 1^1 (x, ,t,a)da = <|)(t)dt 

la 




II 1 { Xi , t ,a)da 
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where e is a positive small number. Using (4.18) for the 
first part of the integral, we obtain: 

a 

^11 (x-j ,t)(j)(t)dt = 

-a 


Making use of the singl e-val uedness condition 

' ® f ^ c 1 1 

(t)(t)dt = 0, the unbounded integral — — da drops out, 

J-a ^0 “ 

leaving only bounded integrals which can be evaluated numer-- 

fb *. 

ically. Using (J>{t)dt = 0 for the second crack, similarly 

-b 

one can show that the singularity due to 1/a cancels In 
all the integrands I... • 



I (j){t)dt I da 

-a ° 

f“ 

+ 0(l)da + I , , { X-. , t , a) da 

J r» J c * 




5. CASE BROKEN LAMINATES 

This is the case when one of the cracks touches the 
interfaces (i.e., a = h^ , or b = h 2 ). The integral equa- 
tions (4.16a,b) are still valid but some of the kernels 
k.jj are no longer bounded. For example for a = h.j k^ 2 » 
*^21’ *^22 bounded but k-| -j becomes unbounded as and 
t approach the ends + h^ simultaneously. In this case 
the integrand I ^ ^ of k-j -j di verges as a 

In order to obtain the proper singularity at the 
crack tips and to compute k^-j numerically, the singular 
part, k^i 5 » should be evaluated in closed form. In this 
case, the kernel k-j^ can be written as: 

kii(Xi»t) - 

where is the singular part and k-j^^ is the bounded 

part of k^ ^ . 

Following the procedure described in [4] ■] g ( ) 

is obtained as follows: 

TTkiisC^i.t) = 

+ A 

86 

+ A 
87 
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(h^-t)/^/|wij + Iwilh^ 

[(■h^-t)v^/Iw,l + |wi|h^]2-(w^x7)^ 

{h^-t)»^/|wi! + Iwalh^ ^ 
[(h^-t)/^/[wiT+ (walh^l^-Cw^x^)' 

(h^-t)/B7/iw3i + |wi|h^ 
[(h^-t)v^/[w3| + !wi|h^]^-(w^x^)' 


i 


+ X 


88 


{[(hi- 


(h]-t)/^/|W 3 [ Iwslh^ 


t)/^/|w3| + tw3|h^r-(w3X,)^ J 


- 5 < h-j 


(5.1) 


The governing singular integral equations become: 


I \ ^ t-V TTk^^^(x.j,t) 

U I 


-h] 
'b 


(j)(t)dt + 


f"i 


kii(Xi»t) - 


(j)(t)dt 


’U ^ (1 -V V ) 

+ ' Ty^E " h^^l^ -h^<x^<h^ 


b *1 


i f" ^ 

IT I t- 
-b 


ill 


fh 


dt + 


^k2i(x2>t)4)(t)dt + J k22(x2>t)((> (t)dt = 

hi -b 

' * * 


Pj(x2) 


2'Yi4Ey 


- b<X2< b 


(5.2a,b) 


Since in the integral equation (5.2b) the only singular 
1 


term is 


t-X2 


, the power of singularity at the end of the 


internal crack in the second layer is still 1/2. But in 
(5.2a) we have further singular contribution from the ker- 
nel l<i ^5 resulting in a power different than 1/2. To find 
this singularity power y, viB will again use the procedure 
described in reference [14]. Throwing all the bounded 
terms to the right hand side, the singular integral equa- 
tion (5.2a) can be written as follows: 


1 (•'i 


ir 




t-x, * ,t) 


({){t)dt = P(x-|) -h-|<x^<hi 


(5.3) 
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where P(x-j) is a bounded function for all values of . 

The unknown function 4 >(t) can be written as (see [14]) 


4>{t) = 


Fit! 




;5.4) 


where F(t) is bounded and Holder-continuous in the inter- 
val |t{<h^, and 0 <Re(Y)<l. 

Define the sectionally holomorphic function: 


Hz) = i 


fh 


ih 


Then , 
^{z) 


1 m dt - i 

t-Z 71 

1 


1 TTY 


f hi 


lilk 


ITTY 


h-j ( t-h-j )'Y( t+h-j )‘^( t-z) 


F(hi) 


dt 


• _] . 1_,+ V (z) 

(2h-| )"''s1miT (z+h,)''^ (2h-| >’'’510117 (z-h,)’’ ° 

(5.5) 


where 




z+h, I 


, Yq < Re(Y) 


C and Yq real constants. 


Using (5.5), equation (5.3) takes the form: 
F(-h^)cot 7 TY F(h^)cot 7 TY X 35 Iw^l F(h-|) 


/3T 
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1 + 1 




^86 
+ -y- 


F(h,) 


(2h-j )^sin7TY 


Nil 1^3 




1 + 1 


“T 


(h^+x^)^ (h^-x^)Y 


^87 ''^3 


F(h.) 


V37 


(2h-^ )^simrY 


IWil lW3 


I vfl 


Y 


1 + ^ 


(h^+x^)^ {h-|-x^)‘^ 


+ -- 


^8 


F(hi) 


T 


(2hi)^sin7TY 




Vi, 


l^. 


Y 


1 + T 


(h^+x^)"^ (hyx^)^ 


= P^Cx^) 
(5.6) 


where because of symmetry F(h.j) = - F(-h.|). 

Multiplying both sides of (5.6) by (h^+x-j)”^ and let- 
ting x-| = -h-| we obtain the following characteristic ^ 
equation: 

IwJ 1 Nii V 


2 cosTY + Kr , - 

2 ' 
Wi 


|wj IW3T 

1 

1 

1 

.. 


Y 


+ X 


|W, 


w. 


87 


j aj. _ 



X J ^ 

“/ST 

r 2 ■ 

W3 


'W ; 



1 


Y 


= 0 


(5.7) 


where X.'s are elastic constants defined in Appendix A 
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This is the same equation found in [7]. Choosing the or- 
thotropic elastic constants close to isotropic constants 
numerically we find the same singularity power computed in 
[ 6 ] and [9]. The characteristic equation (5.7) can be 
solved numerically to find 7 . For practical orthotropic 
materials equation ( 5 . 7 ) has only one root between 0 and 1 . 
To establish the dependence of y on the material constants 
more accurately* a separate study of equation (5.7) is 
needed . 
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6 . 


CASE OF A CRACK CROSSING THE INTERFACE 


To formulate this problem we will start by using the 
crack configuration shown in Figure 3. In this case we 
have again an internal crack in the first layer* but two 
symmetrically located cracks in the second layer. Using 
the symmetry property of f (t), we can write: 

k.2{x.,t)4»*(t)dt = iQ^k-gCx^.t) - = ^*2) 

-b 

and 


f'’ illti ot = 


’ + J 

t-X„ 

-b 2 

0 

t-Xg t+X2 




( 6 . 1 ) 


Therefore we can write the governing singular integral 
equations, by simply changing the limits of the integrals 
from (0,b) to (c,d) in equations (4.16a,b). Thus, we 
obtain: 


i 

TT 


r ^ . 


1 (x^ ,t)4>(t)dt 


-a 


a 

fd 

c 




’ + 1 


t-X2 t+X2 


[k^2(^l»^^ ■ 2^ (t)dt 

* r® 

4) {t)dt + I k2i (x2,t)(|)(t)dt 
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* 

+ J Ck22(x2.t) - k22(x2.-t)]<l> (t)dt = 


(1-v* V* ) 

I . x y. yx .„_ D (x ) 


c<X 2 <d (6.2a,b) 

By letting a = and d = h 2 we obtain the case of a crack 
crossing the interface. As in the previous case for a = h^ 
and d = h 2 all the kernels k^j Lsvi ome unbounded when ,t 
approach the ends + h^ and X 2 »t approach the end h 2 simul- 
taneously. Therefore to study the singular behavior at the 
Interface and to make the kernels numerically integrable, 
the singular parts of the kernels k.j must be separated. 

The kernels can be written as: 


where k.j^{x^,t) is the singular and k^ j ^{x^ , t ) is the 
bounded part, Following the same procedure used in the 
previous section the expressions of k . . are found as 

1 J b 

follows: 


^^1 Is^^l*^^ ” 


DRIGWALPAGEte 

PB POOR QUALXlx* 


85 


+ X 


^ (hi-t)/3T/lwil + Iwilh-j 

U(h,-t)^/lwi| + Iwilh^J"- (w^x^)"j 

{h^-t)v^/|wi I + [Wajh^ 'I 

7) 




+ X 


+ X 




t)/^/.lWj[.+ lw3lh^]^ - (W 3 X 
(h^-t)»^/lw3l + |wi|h^ 
t)/B7/|w3| + |wi |h-j]2 - (w 
{h^-t)/^/|w3l + Jwslh^ 


[(h^-t)/35/|w3i + Iwalh^]^ - (wjX 
-31- 







93 


f (h2-t)*^/|wT| + jwilhi 

\[(h2-t)»^/lwtl + Iwalhi]^ - 


+ \ 


94 


+ X 


I 

I 


* ^ 95 ! 


(hg-^ 

[(t,2-t)/^/lwT| 

(h ,-t)./^/|wT| + |w3|h, 

[(h2-t)v^/|wtl + lW3lh,]"-Iw3X 

(K2-t).^/|wth |wi|h, \ 
[(h2-t)/if/|wt| + l«i|''iP-(«1>'l)V 

(h2-t)/3|/|w 




TTk2is(x2’^^ ’ ^101 


+ X 


102 


+ X 


103 


;y|wt| + Iwajh^ ] 

[{h2-t)/^/lw1|Mw3lhT]^-(w3)^"J 

f (h-|-t)/35/iwi| + lwilh2 1 

\ [(hi-t)v^/lwi| + lwtlh2]^ - (w^X2^^ 

fe 


{[(h, 


(h^-t)v^/lwil + Iw3lh2 

+ lwtlh 2 ]^ - (W 3 X 2 ) 
(hi-t)»^/iw3l + lwtlh2 


+ X 


104 


t)/37/lw3l + lwilh2]"-(wiX2) 

( (h^-t)/3s/lw3l + 

\[(h^-t)/^/|w3l + Iw3lh2]^ - (W 3 X 7 


TTk22s(^2’^^ " ^109 


+ X 


no 


+ X 


111 


+ X 


I 

1 

1 


(hp-t)/g*/lwi| |wi|h2 


[(h2-t)/i*/lwthlwtlh2]^-(WiX2) 
[(h2-t)/^/lwt| + lwt|h2]^ - 


(h2-t)/^/lw3 


[(h2-t)v^/lw3l + jwi 

(h2-t)/ ^/|w*l + lw'llh2 
[(h2-t)/^/l^J| + lwalh2]^ - (W 3 X 2 ) 


+ lwilh2 \ 

;ih2l^ - (WiX2)"j 

■1 
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Separating the singular parts, equations (6.2a,b) take 
form: 

u r* “1 

(j)( t)dt 


the 


and 


fh 

1 -X -1^ 

( X _ , 1 1 

-hi 


lls^’^l * ^ 


+ f kkii(xi.t) - k^T5(x^,tS().(t)dt 

ihi 

+ I ^k*2g(xi ,t)<|)*(t)dt 


+ f'’^[k*2(x,,t) - k*2s(x,.t)]4.*(t)dt 

J r 


'14S 


zy7x !^ -hi<Xi<hT 


i [' 

^ J- 


■1— j — + 7rk««_(Xo,t) 
t-Xo t+X2 22s 2 


(j) (t)dt 


f ^[k2^(x2.t) - k2is(x2.t)]4>(t)dt 


-h 


1 


f ’k2is(x2.t)cj)(t)dt 

-h 


1 


* 1 kjCxj.t) - k*2,(x2.t)]^(t) 



where 




To find the proper singularity power 3 at the inter- 
face, we will first throw all the bounded terms to the 
right hand sides of the equations as it has been done in 
the previous case. Then, we obtain the following system 
of equations: 


1 P 1 

u J tPx7 ^'^1 Is^^l 

-h ' 


1 


4> ( t ) dt 


'2 * 


^12s^^l*^^^ (i')dt - Q-](xi) 


- h^ < < h^ 


IT t-x^ t+*x2 ^’^22s^^2*^^ 


(j) (t)dt 


fh, 

+ j [k2ig(x2»t)(j)(t)dt = Q2(^2^ 


c< X 2 < h 2 


(6.4a,b) 


where Q-|(x^) and Q 2 (x 2 ) are bounded functions of x^ and X 2 

lAr 

Considering the behavior of (f>(t) and tf) (t) at the end 
poi nts , we can wri te : 


<t(t) = 


Fit! 




. ♦ (t) = 


* / \ 
F (t) 


(h2-t)6{t-c)^ 


(6. 5a ,b) 


Define the fol 1 owi ng sect i onal 1 y hoi omorphic functions: 
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( 6 . 6 a,b) 


■Kz) " H"’ ^ ““ 

-h 


hj. .*, 

c 


1 


From [14] and using (6.5a,b) and ( 6 . 6 a,b) we have: 

F(h,) , 


l|j(z) = 


^ (z) = 


F(-h,)e”'B , 


(2h^)^simr3 (z+h^)*^ (2h^)'^simr3 (z-h^) 


^3 


3 o,'. 


+ 'I'oCz) 


.3 "^0 


F*(h 2 ) 


F*(c)e^‘^*^ 1 

(h 2 -c)^simr 6 (z-c)*^ (h 2 -c)^sin 7 r 3 (z-h 2 ) 


1 * / \ 
+ (Z) 


3 


(6.7a, b) 


where 'j^Q(z) and 'I^q(z) are bounded functions which around 


ends behave as follows: 


!z+hJPo 


. 3 o<Re( 3 ) 


and 


z-h2l^o 


. 3(,<Re(3) 




and 




Iz-cl^O 


, 6q< Re{5) 


^0* *^o* ^ 0 * ^ 0 * ^ 0 * ^0 constants. 

Using (6.7a,b) and following the procedure used in 
section 5, equations (6.4a,b) reduce to: 


cotird = 0 


(6.8) 
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and 


F(h^) 

2(2h^ )^simr3 


F*(h 2 ) 

+ — — -r — 
2 (h 2 -c) sinir3 


- 2costt 3 + Xgg 




1 i".i 

^ 86 /B7 I w.l lwsr ^ 

^ m . 

k J V 





+ X, 


|w*| 

■ 1 

+ X l"*l - 

1 


Ps 

|wtl|wt 

6 ^112 ^ 

W3 

3 





= 0 


(6.9a,b) 

Equation (6.8) gives the expected 6 = 1/2 singularity power 
at the crack tip. 


(6.9a, b) is a system of homogeneous linear equations 
for F(h^) and F*{h 2 ). Since F(h,) / 0. F*(h 2 ) t 0. 

3 f 0,1 to solve the system one should equate the determ- 
inant of coefficients to zero. Thus, 


A{3) = 4cos^'it3 - 2costt3<X 


(■ 


* 


1 




w, 




+ X 


w. 


111 




J— + X 
* -3 ni2 


w. 


1 




w, 


+ X 

^85 


1 


/b; (rii) 


'86 


IWi I 


+ X ^ 
8^ ^ 




+ X 


W 3 I 1 


^ (r33) ^ 
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^^iGUslAjj pArp rc 




113 ■; Tnj IH (r r* 1® 

*'^ll''ll' '’^11^3' 


+ X ^ ' 


C-ll^)' ^ <^13^1>' 




1 


+ 

1 



+ 

(rijrii) 

^118 F * vP 

^'^13’^33^ 

('•33’'l1> 

P 



1 


= 0 


(6.10) 


^120 

f * \3 

121 (r r 1^ 

''^33'^33l 



where 

2 

_ ^ 




2 

li- 


’"ll 

' /rr 

’ '^13 ^ 

F8 s 

' 33 

/F7 


if 

w*^ 

1 

* 


* 

*2 

W 3 


^11 


’ 

^ ■ 

33 



From the 

characteri Stic 

equation (6.10) we 

can determine 

the singularity power p. 

Choosing the orthotropic 

elastic 

constants close to Isotropic elastic 

constants, we 

recover 


the singularity power found In [8] and [9]. Equation 
(6.10) does not always have a root between 0 and 1. For 
some orthotropic material combinations, there Is no power 
singularity at the Interface. In this case one should In- 
vestigate the possibility of pure Imaginary or complex 
roots. Numerical computation shows that there are no pure 
imaginary roots or complex roots for which the real part 
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is between 0 and 1 . On the other hand F(h^) and F^Chg) 
are related through ( 6 . 9 a) or ( 6 . 9 b). This is a condition 
to be used while obtaining the solution. The absence of 
power singularity for some orthotropic material combina- 
tions may be very important from the view point of design 
applications. Therefore we will study in some detail the 
behavior of the crack surface displacement derivatives and 
the stresses at the interface. 

Let's first investigate the possibility of a weaker 
i.e., logarithmic singularity in the crack surface dis- 
placement derivatives at the interface. Suppose that the 
power singularity 3 at the interface is zero. Define: 


*(z) = dt . 4>*(z) = i dt 


'-h 


The behavior of 4^(2) around z = + h^ , and of ip (z) around 
z = hg can be expressed as: 

* <l>(hi) 

(z) = — - — log(z-hi) + 'f>o](z) near z = h^ 
4 >(-h,) 

ip (z) — log(z+h|) + <i>o2^z) near z = -h^ 

* 4>*(h2) 

^ (z) = - — - — log(z-h2) + near z = + h2 

(6.11a,b,c) 

where <}>qi(z), bounded functions. 
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Using (6.11a,b,c), near x-j = and equations 

(6.4a,b) take the following form: 


log(h^-x^)l 


<t>(hJ 


IT 


1_X I'llL - X -tei-X - X ^ 

'^85 O /TT— " '^86 O fn — ■%o/^"^88 


w, 


2»W 


2 / 3 : 2/^ 


2/^ 




IT 


W 




w 




93 


2/^ 


94 


w*| |w*i 
3 ' \ 3 > 

95 % r-r Ne 


1/^ 

P 5 


2/^ 

Ps 


- F^(x^) 


log(h2-X2)|~r- 

(j) (h«) 


w 


w 


2/3? 2/37 2/^ 


W I iw I 

^103' ^104 


2/37 


IT 


W 


1 ^ X, 


w 


X. 


w 


w 


:i 

— X. 




F2(x2) 


(6. 12a,b) 


where Fi(x-|) and F 2 (x 2 ) are bounded functions. 

In order that equations (6.12a,b) be bounded for x-j = 
h 2 » the coefficients of the logarithmic terms must be zero. 
Thus: 


4>(hT ) 


w 


w 


!w. 


w 


1 \ 1 1 -\ 3 3 

' ^85 o '^86 o /;r- ^Q7 88 


(})*(h2) 


IT 


2/37 

* . 


2/^ 2/^ 


2/^ 


iwj |wj 

2/^^^^ 

P5 P5 P5 P5 


= 0 


and 
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|w,l N,l |wj |wj 

2./B7 ' 2^ ■ 2^ *103 ■ 2^ ^104 


(j)*{h2) 

+ =- 

IT 

(6.13a,b) 


,1 KL 

2-^ 


w 


w 


'IO'TTtf in 
2 / 3 * 2 /^ 


w 


2*^ 


'112 


= 0 



(6.13a,b) is a system of linear equations for <J){hi) and 

4r iic 

^ (h 2 ). Since and cj> (hg) are different than zero, 

in order to have a solution the determinant of coefficients, 
A, must be zero. Numerical computation shows that A = 0 
and either from equation (6.13a) or (6.13b) we have: 


(h^). 

?{h2) 


and using the symmetry condition <j> (h 2 ) = - <{> (-hg). we 
obtain: 


4>(h^) 


(6,14) 


Relation (6.14) shows that the surface displacement deriv- 
ative is continuous at the interface. This is an important 
result which makes the solution of the singular integral 
equations easier. 


To study the behavior of the stresses, let's first 
write their expressions at the interface. By making use of 
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(3.11) and (4.15) we obtain: 


g(y) ^ g {h^,y) = ^K^(n,y) 0 (n)dn + |^%(n.y)4)*{n)dn 


x{y) = 




xy 


1x ' ihi 

rh 

-h 


frin if 

f— = j K 3 (n,y)<t)(Ti)dn + K^(n,y)<|) (n)dn 

( 6 . 1 5a , b) 


1 


where K.(n»y) (j = 1>4) are given in Appendix B. 

J 

The kernels Kj become unbounded as y -»■ O'*’ and n t 
or n-^-hg respectively. When equation (6.10) has a root, 
i.e. when the functions !j)(t) and i^) (t) are singular, the 
stresses have the same singularity power as we will see in 
the derivation of the stress intensity factors. But it is 
necessary to know whether the stresses have a logarithmic 
singularity when the crack surface displacement derivatives 
are bounded. To do that let's first separate the singular 
parts of kernels K.. Following the usual procedure, we 

vJ 

obtain: 


K,5(n.y) 


4y 


13 


Yi(h^+n) 


(hi-n) 


2 2 
w?y 


- Y 


3s 


(h-i+n) 


1 w?y^ 


tr * 



h,+n _ y,, h,-n 

^ ~~T~2 ~ Xn y, 'i' 2 “ 
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K25(n,y) 


K35(n,y) = 


K4s(n.y) = 


, (h,-n)^/|w 


2 2 


C{h^-Ti)/^/lwil] +y 


1 ^ (h^-n)»^/|w3l 

* ‘^3^83%^84) 


1 , , (h2"Ti)/^/|w* 


2Yi.X„, 3 89 4 90 ^^2 


'13'^SO 

1 

^^1 3^80 


(Y3>^91+Y4^92) 


{h2-n)/^/lw3 


* I i2 . 2 


[(h2-n)/if/|wj]‘ + y 


^11 


wily//07 lwxly/>^ 




T ■*■ 2" 7' 

wiy 


IwalyZ/BT Iwsly//^ 


3s 


” ””1 2 '~ 

2 w y 


+ (h-i+n) + (h^-n) 


' <V81*'^10^82> [(. -n)^/|w.|]S/ 


" ‘^9^83*'fl0^84> 


a n^qn^ ^ . 1 ^ * I t2 ^ 

2^13^80 


^ '"'9’'91 *1^10^92) [(hg-M^/lw’n' + y- 


★,^2 . 2 


(6.16a ,b,c,d) 
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Keeping only the singular terms, equations (6.15a,b) can 
be written as follows: 



where A(y) and B(y) are bounded functions 


Define: 


^(z) = 




-h 


1 Mil dt and 
t-z 

1 


'(z) ■ f 4^ dt 


Considering their behavior around ends, we can write: 


ij;(z) = 4 >(hi)log(z-h^) - 4 )(-h^ )log(z+h-| ) + 4 'q(z) 


and 


i|i*(z) = (j)*(h 2 )log{z-h 2 ) + ^*(z) (6.18a,b) 

where ’I^q(z) and 4 'q(z) are bounded functions. Making use 
of (6.18a,b) and {6.14), equations (6.17a,b) become: 


cf(y) = 


log y 4>{h^) 

1 Iwa 


^ '"2 - TO ^ <^3^81^ V 82) 


1 


iw, 


2^13^80 /F? 




1 iw; 


2 ^ 13^80 ^ 


- (T3^gi+Y4^92) 


+ C(y) 
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and 


T(y) = D(y) (6.19a,b) 

where C(y) and D(y) are bounded functions. Equation 
(6.19b) indicates that the shear stress at the interface 
is bounded. Since C(y) is bounded and (l)(h^) f 0, if a(y) 
is bounded the coefficient of logy in equation (6.19a) 
should be zero. Numerical computation shows that the 
above mentioned coefficient is identically zero. Therefore 
the normal stress a at the interface is also bounded. 

These are important results for orthotropic materials and 
may have practical implications in designing with composite 
materials. 
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THE SOLUTION AND THE RESULTS 

Since we have mainly three different problems, the 
solution will be discussed in three sections. 


7.1 Case of Internal Cracks 

In this case we have to solve the system of singular 
equations (4.16a,b). Defining = aK| , t = ax for 
- a< , t <a and Xg = aK 2 » t = br for -b<X 2 » t< b after 
normalization, equations (4.16a,b) take the form: 


, rl A (t) rl _ fl ^ * 

7 dx + a J k°^(K|,x)<|)Q(x)dx + b |^k° 2 (K| .T)(})Q(x)dx 




^ '1 1 
~ dx + a k 2 i(i< 2 »T)(J)Q(x)dx + b 




★ ★ 


7^P2(<2) 

2^1 4^y 


where the index "o" denotes the normalized quantities. To 
get the complete solution we need also the single-valued- 
ness conditions: 



(7.2a,b) 


Since 4>q(x) and 4>*(t) have a power singularity -1/2 at 
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the ends, the solution will be sought in the foris: 


4>o('f) = 




F„(t) 


— and (frnC'f) = ~ — 

/T^ ° /TT‘ 


where Fq(t) and F^(t) are Holder continuous in the interval 
- 1 < T < 1 . 


Using the method described in [11] we obtain 


N 

J=1 






+ I btlk^2(K,,Tj)F*(Tj) 
0 • 




J ank°,(K,.Tj)Fg(Tj) * 

J I 0 • 


biTkL{ic.,T.) 


T.-K^ - 


* * 

, ■ N ^ . ' . yyx^ . P°(^ ) 1 = 1....N-1 

N „ N „ , 

•I S ° .1 ° 

3”1 3”1 


where 




lb 




Tj = cos ^ (2j-l ) j = 1 , . . .N 
i = 1 , . . .N-1 


K. = cos y 


it 

The 2N unknowns F^j(tj) and F^(Tj) can be found by solving 
equations ( 7 . 3a , b,c ,d) . In this problem we are mostly in- 
terested in the computation of the stress intensity factors. 
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The stress intensity factors may be expressed in terms of 
the density functions Fq(t) and F*(t) as follows: 


For a < h^ : 
arid for b < h 2 : 



lim /2{x, -a) aT„(x,,0) 
xi-^a ' ' 

lim^ i^2{x«-b) cTo„(x«,0) 
X2^.b ^ ^ 


(7.4a,b) 


Making use of equations (4.16a,b) and definitions (7.4a,b), 
after lengthy algebra (see Appendix D) we obtain: 


and 




F„(l) 





F.(l) 


(7.5a,b) 


The computation is done for generalized plane stress case 
only. Results can easily be obtained for plane strain case 
by redefining the elastic material constants. In the per- 
turbation problem considered p-| and P 2 are constant. As- 
suming that there is no constraint in x- direction, p^ and 
Pg satisfy the following condition: 

!i - !/ 

P2 e; 

where and are the Young's moduli in y direction. 

Two material combinations are formed among the follow- 
ing three materials. 
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1. = 5.62 X 10^ N/m^ . = 17.37 x 10^ N/m^ . 

X y 

= 0.492 X 10^ N/m^ . = 0.036 

xy xy 

(E = 8 X 10® psi) , (E = 24.75 x 10® psi) , 

X y 

= 0.7 X 10® psi) . (v^y = 0.036) 

2. E, = 13.71 X 10^ N/m^ , E„ = 3.16 x 10® N/m^ , 

X y 

G^y = 2.46 X 10^ N/m^ , = 0.650 

xy xy 

(E = 19.5 X 10® psi) , (E^ = 4.5 x 10® psi) , 

X y 

(Gyy = 3.5 X 10® psi) , {v^y = 0.650) 

3. E = 15.76 X 10® N/m^ , E„ - 16.78 x 10® N/m^ , 

X y 

G^^ = 6.09 X 10^ N/m^ , = 0.300 

xy xy 

(E = 22.447 X 10® psi) , (E = 22.6 x 10® psi) , 

X y 

(G^y - 8.655 X 10® psi) , = 0.300) 

As it is seen from the values given above the first two 
materials are orthotropic, while the third is isotropic. 
The following pairs of materials are used: 

Combination I: The first layer is of material 1, 

the second of material 2. 

-A9- OEIGINAL PAGE Ife 

OF POOR quality 



Combination II: The first layer is of material 

3, the second of material 2. 

Choosing the same materials and letting a, b, h^ or ^2 go 
to proper limits we recover all the special cases done in 
[5]. [6], and [10]. 

Figures 4-12 show some of the calculated results. In 
Figures 4 and 5 the stress intensity factors k_ are plotted 

a 

versus h 2 /h-| for b = 0 (there is no crack in the second ma- 
terial) and for the two material combinations. For h 2 = 0, 
we recover the solution of colinear cracks imbedded in a 
homogeneous material (see [10]). It is important to note 
that in the colinear crack problem the material doesn't 
have to be isotropic. As h 2 -^“» k^ reaches an asymptotic 
value which can be found in [5]. For a fixed h 2 /h^ ratio, 
kg increases as a/h-j becomes larger. 

Figures 6 and 7 show the stress intensity factors k|^ 
for the case a = 0. In this case also, for h^ = 0, we ob- 
tain the solution of colinear cracks. There is a critical 
value of (h.j/h 2 ) for which the stress intensity factor 
starts to decrease as the ratio b/h 2 increases. For the 
examples done this critical ratio is between 0 and 0.5. 

For h-j the stress intensity factor k|^ reaches an asymp- 
totic value which also can be found in [5]. 
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The stress intensity factors and when both 
layers contain cracks, are given in Figures 8-11. 0 

as (b/hg)-^!, since the power singularity y is less than 
0.5 when the crack in the second material touches the in- 
terface. Another interesting result is obtained from the 
comparison of isotropic and orthotropic materials. As it 
is seen in Figure 12, for the same E and E* the stress in 

J' 

tensity factor k^ for orthotropic materials can be larger 

a 

or smaller than the stress intensity factor k_ for iso- 

d 

tropic materials depending on the other elastic constants. 
One can significantly reduce k_ by a convenient choice of 

u 

the elastic constants. The materials used in the compari- 
son are given in Table 1. The dependence of k^ on the 

d 

★ 

materials constants is given in Table 2. 6 and 6 are 

* xy xy 

the most important constants, while keeping E and E con- 
stant. To reduce the stress intensity factor k_ , it is 

d 

"ft "ft 

sufficient to increase E^, or decrease E^, 

X xy xy X xy 
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7.2 Case of Broken Laminates 


The solution will be obtained by solving equations 
(5.2a,b), with the singl e-val uedness conditions. 



Defining again, 

t = , x-| = h^K^ for -h^<x^,t<h^ 

and 

t = , Xg = b<2 -b<X2»t<b 

the normalized form of equations (5.2a,b) and(7.6a,b) can 
be written as follows: 


1 

IT 


+ TTh ] 


rl 


<{. o ( T ^) dTi + h , J ^ k °^^{ Kl . T ^)<|. Q ( T ^) dT ^ 

0 -V V . v ) 


+ b |^k°2(Ki.t2)<}>*(T2)dT2 = - 


1 '^ o ^'^ 2 ^ 

i ( 0 . ■ dT , + h , 

To-Ko 2 1 


IT 


-1 


L 2 r ^2 


rl rl Q * 

k2i(K2»^l ^‘)^o^'^l^^'^l ^ ^•<22 (k2»'^2^^o^'^2^^'^2 " 


* * . 

-- V V P°(K;) . -1«2<1 


I ' t ’ o *'! 


Ti )dx-j =0 and 


1 


il 




( 7 . 7a ,b ,c ,d) 
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To obtain the solution, we will use the numerical mfttrod 


described in [11]. Hence, we obtain; 


1 N 


■j = l 


-1— + . ,T^ .) 


1 j-'^li 




* ^1 j^’^llf^'^l i ’"^l j^^'o^’^l 

N 

^ ^ j = l'^12^'^li*^2j^’'o^'^2j^ 




1 ^ 




lA'-y 

+ nbk^2<'=21-''2j) 


fo(^2j) 


* ,L "j'^lj>''21<'=2i’''lj>fo(''lj> 

•J • 


* * 

P^(-2i) . 1 = 1.---N-1 

2y,4E^ ^ 


( 7 . 8a , b , c ,d) 


where 


* Yr > - 4*(, 1 - !oil2l 

'Tn^ ’ %<^2> - /p-r 

^ * to 




= 0 , J = 1....N 

= 0 , 1 = 1....N-1 
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"^2j = cos ^ (2j-l ) , j = 1 N 

7T i • 

**'2 i ~ cos -jq~ » i = 

and Wj(x^j) are the weights of j ^ • 


Solving the 2N x 2N system of linear equations one 
can find the 2N unknowns F^CT-jj) and f^o(T 2 j)- But again 
we are interested in the stress intensity factors, which 

•fg 

can be calcul ted in terms of and 


Define: 


and 






02y(x2 ,0) 


After some calculation shown in Appendix D, we have: 






(7.9a,b) 




★ . . 

fo<’) 


The results for the case of broken laminates are shown in 
Figures 13-16. Again the same material combinations are 
used. In Figure 13 the stress intensity factor k, is 

a 

plotted versus the ratio h 2 /h-j for b = 0. When h 2 -*■ k^ 
has an asymptotic value which can be recovered in [7]. 
Figure 14 shows the variation of with h-|/h 2 , for the 
case a = 0. For b 0, the variation of k^ and kj^ with 
b/h 2 are given in Figures 15 and 16. 

7 . 3 Case of ^ Crack Crossing the Interface 


In this case the governing singular integral equa- 
tions are (6.3a,b). As it was pointed out in Section 6, 
the characteristic equation (6.10) does not always give a 
singularity power at the bimaterial interface. Therefore, 
the numerical solution needs care, and we should solve 
equations (6.3a,b) considering the singular and non- 
singular cases at the bimaterial interfaces. 

7.3.1 Singular Behavior at the Interface 

For the material combination II (isotropic-ortho- 
tropic) equation (6.10) has a root between 0 and 1 . Using 
Newton-Raphson method to solve equation (6.10), we have: 

B = 0.04248 . 
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We win make use of the following definitions to normalize 
equations ( 6 . 3 a,b), the s ingl e-val uedness condition, and 
the relation ( 6 . 9 a) : 


t = h^T-j , = h-jK^ for -h-|<t, x^<h^ 

h«-c h«+c h«-c h2+c 

t = Tg + \ = “7~- <2 


Then, we obtain: 


-1 ^ ^ -1 


( l - V ^„ Vyy ) o . 


M ’ 

"'-1 




2 '" 2 ’ 


♦„( T 2 ) dt 2 


h, |’k2is(<2’''l>*oh>''^l * h|,'‘21f(''2>’'l>*o‘^l>‘‘''l 

h^-c fl 


-1 -1 


P^(<2> -K''2<1 

2Y-|4Ey 


/ ^ 


(|)o(Ti)dTi = 0 and F*(l) = - 


il 


hi 


hg-c 


B a 


( 7 . 10 a,b,c,d) 
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where 


a-j = - 2cosTr3 + Xgg 


jwj 1 . , L 


+ X 


v^ 3 T 

1 — r' 
Wi 

S ^86 /PT 


a. 


IW3I 


1 * li 




< ) 


W 3 I 1 


87 




twjk 


vC 


^ 88 


w, 


“2 ' ^93 


Nil 1 


W.ll"l 


. ^ 


. 1 ^ 
B 9 -* /gf 


1 




►IT 


+ X 


95 


1 w* 1 W 3 I 1 

i '^3 ! + X — — 








B "96 /g*- 


w.ll< 




Using the numerical method given in [111 equations ( 7.10 
a,b.c,d) further reduce to: 


N 


n . 1 , 

J = 1 


+ I Wj(t2^) 

j=T 




-I— + TTh^k°Ts(Kii.Tij) + 






h?-c ,*o , V 

2- -las^'^li'^Zj' ~T' "l2f'''vi’''2j' 


k^L(i<n-X;n) + 


*/ \ 


<^;V aIp°(K,,) i = i....N-i 
- 57 - 


1 N 

w I 


j=l 




'2j 


‘‘^21 ^(hp+c)^^ T~ ‘^22s^'^2r'^2j^ 

To 4 +K^.-+ -ir~ 


^23 “R^c 


^2“^ *0 


+ TT -j— k22f(K2i*'^2j^ 


Fo('' 2 j) + I «j(''lj)['’l'‘? 2 s('' 21 -"lj) 

iD * 


M * * \ 

+ Rik2if(K2i»'^lj)]Fo^'^lj^ = ~~ T^r P2(i<2l) i = 1.--.N-1 

N 


and 


where 


'’I J, = ° 


r ^ 


F^O) = - v7 


hg-c 
i. j 


32 0 


F„(l) 


* (.) . iollL , **(,) . 

“ (i-t^)B O-Tjeo+T)*^ 


(7.na,b,c,d) 


P;-®'-®)(T, .) = 0 


j = 1,...N 
i = 1....N-1 
j = 1 ,...N 

Pn-/’'®’‘‘^^ 21> =0 i = l,...N-l 

corresponding weights of 


•2j 
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Solving the 2N x 2N system of linear equations, we 
obtain the 2N unknowns stress in- 

tensity factors can be defined as follows: 



By making use of definitions (7.12a,b,c) and after lengthy 


calculation shown in Appendix D we obtain: 


1-ft * * 

2' 


If = . 'T / 

° (1-M V ^ 

' xyyx' 


f;(-i) 




®xy . J 

f-i® 






'n 






IWal 


/3s 


(Yg^si+Yio^az^ IwJ ^Y9^83'*’^10^84 


o 

00 

kJ 

o 

00 

r< 



jWli 

kJ 



Fo(-l) 


(hg-c) 


^ 13^80 


(Y9A89+Yio^9q) 

Wi! . ^^9^91 ■^^10^92^ 

]wIL 

F^O)- 


r 1 * 1 1 

KL 
M . 

3 

•-i? 

[Jill 

[*'35 J 

3 



( 7 . 1 3a ,b ,c) 

Extrapolating the results found from equations (7.11a,b,c, 
d) the stress intensity factors k|^, k^^ can be com- 

puted in a straight-forward manner. The results are shown 
in Figures 17-19. Figure 17 shows the variation of 
with c/h 2 » for different values of (h-|/h 2 ) ratio. k|^ in- 
creases as (h.j/h 2 ) increases. Figures 18 and 19 show the 
variation of k and k with respect to c/ho. 

A A A ^ ^ 

7.3.2 Non-Si ngul ar Behavior at the Interface 

In this case, the characteristic equation (6.10) has 
no root and therefore the surface displacement derivatives 
are bounded. Since, as it was shown in Section 6, the dis- 
placement derivatives are continuous at the bimaterial in- 
terfaces, using the singl e-val uedness condition to write 
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the integrals from 0 to il, equations (6.2a,b) take the form 
(see [12]): 




r-s r+s 


6{r)dr + k(r,s)6(r)dr = p(s) 0<s<Jl (7 

J n 


where 


with 


G(r) =< 


(}){t) (0<t<hp 0<r<h^) 

(-h 2 <t<-c, h^<r<Jl) 


0<s<P,) 

P(s) * * 

* * (■■h 2 ^X 2 <“C> h-j<s<il) 

2Yi4Ey 


k^{(x-],t) (0<Xp t<h^ , 0<r, s<h-| ) 

ki2(xi ,t) - 1 


k(r,s) 




' + ’ 


k* (x M - i 

X2i\X2»u; ^ 


t-X^ t+X-j 


' . ' 


0<x^<hp 0<s<h^ 

-h2<t<-c, h-|<r<5- 


t-X2 t+X2 


-h2<X2<-c,h^<s<il 
0<t<h-| , 0<r<h^ 


k 22 (x 2 *t) (-h 2 <X 25 t<-c, h-|<r, s<«,) 


k^.j(x^,t) = k^j{x^,t) - k^.j(X|,-t) (i.j = 1,2] 


ij 
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Now we have the governing equations for a crack imbedded 
in a non-homogeneous material, obviously with power singu 
larity -h at the crack tip. Normalizing equation (7.10) 
by means of: 


r = At and s = Ak 


we obtain: 


M’ 

” Ja 


- + 4- 

T+K 


T-K 


1 

G (j)dT + a| k (t,k)G (x)dT = p^(k) 0<k<1 (7.11) 


Equation (7.11) reduces to a set of linear equations 
by using the method of collocations (see [11]): 


N 

I 

j=l 




Ho(Tj.) = 2NPq(k^.) i = 1....N 

(7.12) 


where 


T . 

J 


= cos 


23-1 


•m 


K. = COS 


[A 


S„(T) 


Ho(t) 


/TT72- 


3 = 


. 1 = 


1 , . . . N 
1 , . . . N 


The N unknowns found from equation (7.12) in 

a Straight-forward manner. 


Defining the stress intensity factor at the crack 
tip as : 
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'‘b ° b2y(*2.0) 

we obtain: 

★ ★ _ 

„b . ■ H„(1) (7.13) 

''’xy'’yx 

Using the same material for both strips, we recover again 
the results of colinear cracks in homogeneous medium. Fig- 
ure 20 shows the variation of kj^ with c/h 2 * k^^ increases 

as h^/hg increases. 
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8 . 


CONCLUSIONS 


The fracture problem of layered orthotropic composite 
plates has been studied. The following results have been 
obtai ned : 

1) Depending on the elastic constants, orthotropic 

materials can be classified in two groups: materials of 

type I and materials of type II. (A different formulation 
is needed for each combination.) 

2) The colinear crack solution is the same for homo- 
geneous isotropic and orthotropic materials. 

3) In the case of an internal crack in the first 
layer, the stress intensity factor k^ can be reduced sig- 
nificantly by a proper selection of the elastic constants, 

4) For the case of broken laminates there is a sing- 
ularity power which can be found from equation (5.7). The 
singularity power y varies between 0 and 1 for different 
material combinations. 

5) For a erack crossing the interface, the singular 
behavior at the interface disappears for some material com 
binations. In this case the crack surface displacement 
derivatives are bounded and continuous, and all stresses 
are bounded at the bimaterial interfaces. 
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9. RECOMMENDATIONS 

In the present work, a general formulation of the 
fracture problem of layered orthotropic composites with 
periodic cracks is given. The formulation is done only 
for the case where both materials are of type I. Follow- 
ing the same procedure, the problem can also be studied 
for orthotropic materials of type II, or for the combina- 
tion of type I and type II. The dependence of the singular 
behavior at the interface on the elastic constants can also 
be investigated. 

In our formulation the thickness of the adhesive 
bonding the layers has been neglected. The study of the 
effect of the adhesive also can be recommended. 

A more realistic approach also would be to study the 
problem of finite number of strips. But this problem re- 
quires lengthy algebra. 

There are many other problems to be studied in the 
fracture of composites. We hope that our work will have 
a small contribution in the study of these problems. 
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Table 1. 

The elastic constants 
and material 
combinations 

I 

"j 


Mat No. 


1 

5.62 X 10^ N/m^ 
(8 X 10^ psi) 

2 

13.71 X 10^ N/m^ 
(19.5 X 10^ psi) 

3 

15.76 X 10^ N/m^ 
(22.447x10^ psi) 

4 

17.34 X 10^ N/m^ 
(24.3 X 10^ psi) 

5 

1.025 X 10^ N/m^ 
(1.46 X 10^ psi) 

6 

3.09 X 10^ N/m^ 
(4.4 X 10^ psi) 

Comb. 

No. 

Material of 
the first strip 

I 

1 


II 


3 






























Material 

constant 



Ex 

increases 

increases 


i ncreases 

decreases 


i ncreases 

increases 

IBSlHi 

i ncreases 

decreases 

^xy 

increases 

increases 

V* 

xy 

increases 

decreases 

E and E* are kept constant 

- JL „ , I — — 


Table 2. Dependence of on the 
elastic constants 
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FIGURES 


* 

1-20 








Figure 3. Initial geometry for the case of a 
crack crossing the interface. 





I I I I 1 

O 1 h,/h2 2 


Figure 6. The st»'ess intensity factor for the 
crack in buffer strip (Combination I )• 

- 75 - 


ORIGINAL PAGE fo 
OP POOR QUALITY 











Figure 9. The stress intensity factor for the case 
in which both strips contain cracks 
(Combination I)- 
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Figure 14. The stress intensity factor when the second 
laminate is broken (Combination I and II). 







b/h2 1.0 

or kjj when the first 
e second contains a 
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Figure 18. The stress intensity factor for a crack 
crossing the interface. 
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Figure 20. The stress inters 
crossing the inte: 
at the interface). 
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APPENDIX A 


Definitions of the material constants; 

A superscript * will be used for the material in the second 
strip. The constants 3^, (i = 1,..,10) and Wj (j = 1,..,4) 
are given by equations (2.15), (3.5), and (3.8). 


V w -,39 

Y ,2 - - sign(W 3 ) 6 ,o 

/3s 


Y ,3 = sign(w,) 8 g - ^slgnXwgjBio 

’fa ' “1 ^ ^x »7 

^ - 1 r^5 Ye ^ui 

^14 Y Y 

*13 '13 '12 

* 

^4 ' «3 ^ '’yxSs 

E 1 - V V 

XV yx 

- p * * 

1 1 - ''xy^x 

, ^ ®9“1 

= ''xy ^ ^ 


y,=v 

« xy ^ 

^3 " ^8^3 “ ^7'^4 


* ★ 

^7 ” "’xy"! ^ h 

^4 " ^ 1^337 “ ^7^3 

'fS ° '’xy“3 + ®8 

^5 " ^1^4^7 ■ ^8^3 

Y9=-1 +37Wi 

^6 "" ^9 ■ ^10 

Yio=-’'"V3 

* 

^7 "" ^9 " ^9^2 

Yii = - * "' + sign(w^)3g 

* 

^8 ^9 ’ ^10^2 
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^24 ^10^16 “ Vi 6 


^10 ^ ^3^8 

X^l = X3Y9 

X-|2 == ^6^4 
^13 " ^5^6 
^14 " V3 
^15 " V? 

X^6 = ^337 
^17 ^3^7 '*’ Va 

X^8 ^4^7 

•*f 

^19 ^ ^3^8 '*’ ^5^8 
^20 " ^5^7 
V ^ ^3 ‘ ^3% 
X22 = Y337 
^23 ■ ^7^8 


^25 

" ^10^17 


^26 

= X^ 2 ^^g 

■^10^18 

hi 

= XgX 2 Q - 

■ V 16 

hs 

= Vl9 


hs 

= X^2^19 

“^3^17 

^30 

"^13^18 ■ 

■X^2V 

^31 

" ^11^16 

" ^9^16 


\ - X X 

^32 ■ 11 17 

V ' ^16^12 ’ ^11^18 

^34 ^ V 22 “ ^14^16 
^35 " ^9^21 

^36 ' " ^ 14^17 " ^ 21^12 

hi " ^14^18 ‘ ^22^12 " 
2 

^38 ^ ^15^16 ■ ^7^9 
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Sg V23 


^54 " ‘ ^10^23 

^40 ^5^7 " 

^2^23 

^55 ^13^23 "^5^19 

2 

^41 ^ ^7^12 “ ^ 

15^18 "" ° 

Xgg = - X3X20 

^42 "" ^3^16 


\j ^3^19 

^3 " ^3 ^17 


^58 " ^4^16 ^^8 ■ ^9^ 

"" h^^Q 


^59 ” “ ^5^16 ^’'^9 ' ^7^ 

hs " ^10^16 " 

^11^16 

^60 "^9 ^^4^9 " ^5^7 

^46 " ^11^20 “ 

^1 3^1 6 

^61 " ^3^7 ^^10 “ ^9^ 

\j ^ ‘ ^11^19 


^62 ^ ^8 ^^3^17 ■‘‘ ^4^21^ 
★ - 

^8 ^ ^14^16 “ 

^10^22 

^63= V 

^49 ^13^22 " 

X14X20 = 0 

^64 " ■ ^7 ^^0 ‘ ^9^ 

^50 ^ - ^10^21 


^65 " ^8 ^V23 ■ ^N7^ 

^51 " ^14^19 

^13^21 

^66 " ^ ^^7^9 " ^5^23^ 

^52 " ^10^7 " 

^15^16 

^67 " " ^4^16 


2 


^53 ^ ^15^20 ■ 

■ ^367 = 0 

^^68 5 16 
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^69 "" Vl7 ■ ^4^19 

^70 ^ ^3^ ^7 ^^10“ '*’ ^8^4 ^^8 "^9^^ 

^71 ~ ^3^^8^^n"^9^ ■‘‘ ^8^5^^9 "^7^^ 

* * 

^72 = ^7'^n^5"^3^13) ■’■ ^8^ ^3^1 2 "^11^4^ 

^73 " ^21^ho"^9^ 

^74 ^ ^12^3 “ ^10^4 '’’ ^7^ ^14^3" 

^75 " ^9^5 " ^13^3 '*' ^8 ^^9^3" ^14 ^3^ 

^76 ^ ^23^ ^10* ^9^ 
hi ^ ^7^^0^"^15^3^ 

^78 ^ ^8^ ^15^3“ ^^9^ 

^79 ^ ^3^H^8"^5^7^ 

^80 “ ^24 ^25 ■*■ ^26 hi ~ ^28 ^29 

^81 “ " ^70 " ^71 ■ hi ‘‘‘ ^ 9(^73 ^74 ■'■ ^75^ 

T^ ^^76^ hf h%^ " “Jw^ ^^43"^57'*‘^79^ 
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'82 


” ^58 ’ ^59 ” ^60 ‘’’ ^62 ■*■ ^63^ 


i/ft~Yi y *^s 

+ (Xg^+ Xgg+Xgg) - I ^^67'^W^69^ 


w. 


'83 


^^70'^^71'‘'^72^ " 7^ sign(w2)e^Q(X73+Xy^+X7g) 
12 ' 12 




Y 

1 1 5 


Y w 

12 3 


(X 76 +X 77 +X 78 ) + --y^-]- (^43"^57'*‘^79^ 


^84 "" Y^^ ^^58'^^59’^^60^ ' Y^^ 
12 12 


Y Y Y 

^^64‘^^65'^^66^ ■*■ ■■ji^^^^67'^W^69^ 

la's' ' ’ ' 


Y 7 X 


7''81 


'85 


^1 3^14^80 

^8^82 
^13^14^80 

^7^83 
^1 3^14^80 

^ = ^8^84 

88 Yi3Yi4XgQ 


'86 


'87 


Xon “ — ^ 


* . „* 


^89 


70 


'71 - ^2 - sig"(«i)^g(^3 " ^4 ^ ’'75> 


■ ^1*^ 7^ ^^76^^77^^78> ■ ^2 “|w7 *^^3'^57^^79> 
'1 '1 


* * 


'90 


" ^58 ■ ^59 ” ^60 " )3g(Xg.j + Xg2 + Xgg) 


" j^*j ^^64‘^^65'^^66^ " ^2 J^*j^ ^^67’^W^69^ 
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'91 


*. * 


Y Y 

# (’^ 70 *^ 71 *^ 72 ) * ^ s< 9 "(« 3 >® 10 (’‘ 73 *V 4 *V 5 ) 


12 


4p 

•y ” *Y” 


*•^^1 (V6*’^77*’'78> * '- ‘ | w* r ^* <^43-V’^79> 


12 3 


'92 


12 


Y 


*. * 


1 2 


+ ^ ^^ 64 ‘^^ 65 '^^ 66 ^ '*' ^^ 67 ‘^^ 68 ‘^^ 69 ^ 


12 '"3 


'93 


J^ 7^89 

^ 13 ^ 14^80 

^ 8^90 


94 Yi3Yi4Xqo , 


'95 


'96 


7 7^91 

^? 3 ^ 14^80 

^ 8^92 
^ 3^1 4 ^80 


V; ” ’ ^45 " ^46 ■ ^47 '*’ (>^43 ‘‘' ^50 ’’' ^51^ 


Y, I ■ ^' o 

* -nsT <^52*^54*^55) - i^r <H 2 *^ 56 *' 57 ) 


fw 


^98 ■ ^31 " ^32 " ^33 ‘‘‘ ^^9n{w^ ) 39(^34 + 

Y Y /T 

7^ ^^ 38 '^^ 39 ‘^^ 40 ^ ■ 1" ^ " ^42 ‘ ^43 ^ 44 ^ 


^99 Y^ -Hs ■‘' ^46 ■*■ ^47^ " V~^ ^’9n(w3)3^Q{X^g + XgQ + Xg-j) 
12 12 

' ^^ 52 ^^ 54 '^^ 55 ^ 'iwrr ^^42 '^^ 56 '^^ 57 ^ 
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'100 




(X3^+X32+X33) - sign(w3)p^Q(X3^+X3g+X3g) 

12 12 

Y Y Y 

<’-38'^^39'^^4o| <• ^42-^43'^V4> 


12 ' 3 

* 


^97^7 
^13^*4^80 

^98^8 
^13^T4^80 

^99^7 
^13^*4^80 

= ^100^8 

^1 3^74^80. 


'101 


'102 


'103 


*. * 


'105 


’ ^45 " ^46 ■ ^47 " (X^g + Xgg + Xg^) 


Y, 


v*./^ 


'106 


" jw*l ^^52'^^54‘^^55^ ' ^2 ^\2'^^56‘^^57^ 

“ ^31 " ^32 ■ ^33"Sign(w^)gg (^34’’'^35'^^36^ 

"Y^ Y^ 

“ ^1*^ rir ^^38'^^39‘^^40^ " ^2 7^ *^^'^42‘^43'^^44^ 


^107 "" ^ ^^45'*’^46‘’‘H7^ ■'■ Y*^ ^^9'^^^3^^10^^48'’’^50'‘‘^51 ^ 
12 ' 12 


Y' 


*. * 


Y?7->Sf 


'108 


"ft "ic 

+*^^1 7^7+ ^^52'^^54‘^^55^ "T~*f ^ ^^42‘^^56'*'^57^ 

12 I'^si l^gi 

Y * Y* 

7^ (X31+X32+X33) + sign(w3)a^Q(X3^+X3g+X3g) 

12 12 


Y*^ Y^ 


Y" 


+ Xi /«* ..^ . (X-3Q+X^Q+X^f^) + ^ ^”^42~^43^^44^ 


1 ^5 Y* I *i ^^38 ^39 ^40^ 
12 IwJ 


* 

W 3 
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^109 

^Up 

^in 

^112 

^113 

^114 


^115 

^116 

^117 


ot ^ 


^ 105^7 

^ 13 '^T 4^80 

^ 106^8 

^ 13 ^ 14^80 


^ 107^7 

^* 3 ^ 14^80 

•k 

Mo 8^8 

^ 13 '*’l 4^80 

P 5 

i«ii J«IL X X + 

'' 85^10 

.i^liiiiLx X 


W, 


1 ^ 


^ 93^101 


1^1 1 hiL X X 


w, 


Iw. 


Xi^ f Xi 


IM 

i<l 

X X - I'^^i 

★ 

IM 



W112 ^ 


IM 

ILll 

XX 

i“ti 

/37 


^^86 109 ^ 





■'95^102 


J^lwiix X 


W, 


W 


/z 


^ 87^109 


iL X X 
^ 94 n 01 


w 


i_ 


w. 


^5 

Iw. I [w*| 


^ 87^111 


I^llil^X X +Jliil^x X 

WllO^ r- 5 - Will 

4. l)y.al l^.iJ. X Xi -, r 

.mj^X, 3 X ,04 

.1M^X35X,33 

/87 


Li. 


^5 
, * , 


IW. 




^ 94^102 


WlOl 
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Wll 2 


iw . 


|W 


ni 8 




^ 87^112 


| w . 


iil 


| W , 


iM 


^ 95^104 


| w . 


IW. 


w. 


'119 


lw*l 

^ 88^109 


/ F 7 

I * I 

iM 


Vl 02 


w , 


Wl 03 




|W, 


I * r 

1 ^ 


^ 94^104 


IW 


1. 


IW 


3 j _ \ \ 

^ 96^03 


'121 


/ 3 F m 

X X - i^-i^ X X 
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APPENDIX B 


Expressions of the functions used in equations {4.14} and 
(4.15): 

ra 


1 


Rt (a) * “ 77 

1 3« j 




1 


2^13“ -b 


-(hi-n)a/67/ |wi j ■(h^-n)g»^/ |w3 ! ' 

■ y■^2 •■PvTildn 


,-(h2-Ti)a»^/|wt|_ !lle“(b2-n)av^/!w3ii/(,^)dn 

Yi2 


R2(a) 


1 






- sign(w3)$,.e-<'’r''>“'^/l*'3l 

^12 


<{>(n)dn 


2^1 3°^ 




-:-y-sign(«*3)B:.e-<Vn)“^/l«ll 

^12 


w. 


<j)*(n)dn 




JK- 

^ia 


Xi -(h-|-n)a/^/!wi 




. ^11 ,-(hi-n)g>^/|w,| 


\zW 

rb 




2 Yi 30 t ib 


<j)(n)dn 




iw:i 


(j)*(n)dn 
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R4(“) 


^13“ i. 




1 g-(h-|-n)av' 35 /lWj’ 


w. 


1 g-(hi-ri)a/ 3 r/lw 3 


IW 3 


4 >(n)dn 




2Y*30t L 


Iw, 


I_ p-(h 2 -n)a>^/lwi 


w, 


-(h 2 -n)a/^/|„* 


4>*(Ti)dn 


f(a) = ^24^1 ^26^3^^^ ^27^4^“^ 

•ff •fg 

f^(a) = tanh(w^ah2)tanh(w20th2) 
f 2 (a) = tanh{w 2 ah 2 )tanh{w.|ahi) 

^ 3 ( 01 ) = tanh(W 2 ah 2 )tanh(w 2 “hi ) 

4r 

f^(a) = tanh{Wiah2)tanh(w2ah-j ) 
fg(a) = tanh{w*ah 2 )tanh(w^ah^ ) 
fg(a) = tanh(w^ah^ )tanh{w2cth^ ) 

g(a) = X2-jtanh(w*ah2) + X32tanh(WiOth^ ) + X23tanh(w20ih^ ) 

h(ot) = Xg^f^la) + >^ 35 'f 5 (a) + ^ 36 ^ 6 ^^^ 

m(a) ~ X2gf^(ct) + ^39^5^.®^) ^ 40 ^ 6 ^^^ 

n(a) = - X^ 2 ^®'^^(w-jah 2 ) - X^ 3 tanh{w-|ah^ ) + X^^tanhCWgah^ ) 
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g^Ca) = X^gtanh(w2Qih2) + X^gtanh(w2ahi ) + X^^tanh(w.|ah^ ) 

«io<“) = ^52*^3<“> ^54*’2(“> ^55^6<“> 

no(a) = + Xggtanh(w2ah^ ) + Xgytanh(w^ah^ ) 

ilr if 

g^(a) = Xggtanh(w20th2) + Xggtanh(w-|ah2) + XgQtanh(w^ah^ ) 


h-j (ot) ” ^ ^ 62 ^ 2 ^^^ ^ ^ 63 ^ 5 ^^^ 

nil(a) = '*’. '*’ ^66^5^°^^ 

if if 

n^(a) = Xg^tanh{wgah 2 ) + Xggtanh(w-jah 2 ) + Xggtanh(w^ah^ ) 

it if 

92 (a) = X^gtanh{wgah 2 ) + Xy^tanh(w^ah 2 ) + X^ 2 ^S'^^’(w 2 ah^ ) 
m2(a) = X^gf^(a) + + ^ 78 ^ 4 ^“^ 

if it 

n2(a) = X^2'^®nh(wgah2) - XgytanhCw^ahg) + X7gtanh(w2ah-| ) 


Expressions of the functions k.(j=l,8) used in Eqs. (4.17) 

sJ 


■ 2y^3f(a) 


cosh{wi ax-j ) 

cosh(w,c ' h,) ^7 -92'“> * ) 69 '’ 2 t«‘) 

'I A 


+ " jw 1 ^ ” 2 ^“^ 


cosh{w3ax.| ) 

^ cosiT(w^aH^ ^8 


+ sign(w^) ggh^Ca) + /Bg " “j 




-Qi (a) 

t *r 


w. 


n^(a) 
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1 

cosh{w^ax^ ) 

2Y^3f(a) 

cosh(w^ah^ ) 


Y Y 

~ sign(w3)3^Qli2(a) 
'12 ''12 


Y Y Y 


Y VT '2' 

12 3 ' 


7 ^ sign(w3)S,gh,(a) - ,^^^m,(a) + ^+j-^n,(a) 

12 * 1 o ' a ' a 


Y Y 
11 2 

~T^ 

12 ' 3 ' 


cosh(w3ax^ ) 

^ coshlw^ah-j ) 

Y /r 


1 1 


12 


g^(a) 


1 

cosh(w^ax-j ) 

2Y*3f(a) 

cosh(w^ah-| ) 


*. * 


- 92 (a) - sign(w^)6gh2(a) 


Y* Y* 

■ "'2'“’ ■ ^2 


cosh(w^ax-, ) 

+ ^ - r Y 

coshlw^othi) 8 ' 


*. * 


l 


Y Y ^^^3* 

-9^ (a) - sign(w^)3gh^(a) - - A 2 (a) 

® I w I w j 


1 


^ 4 (^ 1 ’^) " 


1 


2Yi3f(a) 


cosh(w-]aX] ) 
cosh(w.ahTy ’'^7 


Y* Y* 

+ /^ A 


3 c 

® Y_ W 3 ^ 


«1 (j.ii'l 

>* 

u 

! W 


* 


Y* Y* 

S2<“> -r ^'s"<“3Ko'’2<“’ 

, 12 ’12 


Y* A /r* 

(a) + -Ai^^n 2 (a) 


1 2 


cosh(w3ax-j ) 

^ coshC w^ahi ) "'^8 ' 


1 1 


Y* * * Y* Y* 


" in, (a) 


12, 

y* A 


W. 


12 


" 2Y] 3f(a) 


cosh(WaaX2) * 
Y- 


Y, 


f /r -r-S-niJa) - -|^n^(a) 


5 IwJ 0 


cosh(wfahp) ^ 

Y. ./F 
1 


-9o(a) + sign(w^)39hQ(a) 


n 

cosh(w3aX2) * 

•f _ Yg 

cosh(w3ah2) 


-9(a) 
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Y Y 

+ sign(w^)3gh(a) + m(a) - — n(a) 


1 

cosh(w^ax2) 

2Y^3f(a) 

cosh(w^ah2) 




Y Y Y ^ 

■ iwT '"0*“’ "0^“^ 


Y Y 

g (a) _LL . _LL 
^ 0 ' ^ V Y 

12 12 


sign(w3)3^QhQ(a) 


I 


12 


- s 1 gn(w 3 )e,(,h(a) - /T m(a) + n(a) 


cosh(w3ax2) * 
+ ^ Yg 

coshCw^ahg) 


g(a) 


12 




1 

* 

cosh(w^ax2) 

2Y*3f(a') 

cosh(w*ah2) 


*. * 


Y’' . y* ^ 

. 'll P 5 


’ ' Iwtl “ ^ l<l “ 

★ 

* * . . __ Y 


-9o(a) - sign(w-|)eghQ(a) 

-g(ot) 


cosh(w3aX2) *' 

+ Yg 


cosh(w3ah2) 


- sign(w^)0gh(a) - m(a) - X 2 n(a) 

1 Wj 1 Iwj I 

kg(x2,a) = 


1 


2y f(a) 
13 


Q0sh(w-,ax2) * 
T Yy 

cosh(W|ah2) 


9o(“) 


* 


+ s1gn(W3)e,3 h^(a) + X, 

^12 


Y*. Y* 


1 1 


Y*, Iwtl ° 


m^(a) 


Y* X 

+ nja) 


cosh(w3ax«) * 
+ Yc 


lw*l J cosh(w3ah2) ® 

•y'k 'Y*'^ Y 

etnMo) + X, -Ji — in(a) + 


Y Y * 

g(a) + signCwg) < 

y*2 y 


y* X 

u.^-^ n(a) 


12 


The kernels K.(j = l,4) used in Eqs. (6.15a,b): 

V 


.-r ^ 
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cosah^sinanda 


1 r~ 

K,(n,y) = 2v” 

-{ha^ri)a/35/lwil 


^^e"lwJay//07 _ Lj. e"|w3|ay/>^ 


'2 Y 


12 


• 5 T 


fel. 


Y 392 (a) ■•* Y 49 -|(a) 


+ s 


ign(w^) 3 g 


Y 3 h 2 (a) 


+ Y 4 h^(a) 





Y C 

' 

Y3ni2(a) ■♦■Y4nii(a) 

k j 

11 5 

l«.l 

Y3n2(a) +Y4ni(oi) 

^ . 


g-(hi-n)a/37/lw3l 

Y 

f 

" 

Y392(“)'^T49i(“> 

' 

- “ sign(W3)S,(, 


12 

1 

J 

12 

V 


Y 3 h 2 (a) 


74^1(01) 


"Y,, y. 

12 ' 3 


Y 3 n 2 (a)+ Y 4 n^(ot) 


Y 3 m 2 (a) + Y 4 m^(a) 


1 Y ^ 

11 5 


cos 


da 


r 



* 

■'0 

^ V 


K 2 (n.y) = yir 


Y392(^) 


L-(hz-n)a/^/|w*|_ g"(h2+n)a/^/|wi 


*. * 


- sign(w^)3g 


Y3h2(a) + Y4h-|(ct) 


* 

Y ^ 

- 


Y 3 lTi 2 {a) 


Y* 

II P5 i 

i«:i 


+ Y 4 ni| (a) 

-(h 2 +n)a/^/|w* 


Y 3 n 2 (ot) +Y 4 n^(a) 


V+< 


,-(h 2 -n)oi/^/|w* 


- e 


Y 


1 1 


5 

* . * 



* f 

- < 

'Yii 



/ 

' V 


Y 392 («) + Y 49 T(a) 


sign(w 3 ) 3 ^Q 


Y 3 h 2 (a) 


12 
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+ Y4hi(a) 


L 

Y^n^Cot) I 


* *\ 

Y Y ^ 

iwTi 


Y3m2^®^ ■'■ Vi 


I 


Y* 

+ .l.n ^s -2. 


w, 


Y3H2(“) 


TO 


K3(niy) = ^ 


13 


^da 

•00 

« 

-lwilay//3T -Iwaloty/ZeT 

0 "" 

0 111 - e 


/ 


sinah^sinanda 


i_r 

Yi-J Jo 


(hi-n)ct/e7/lWi 


'^13 

y^Qtanh(w3ah^ )g-^M 


YQtanh(Wi0th^)g2(«) 


+ sign(w-j)6g 


+ 

l^ 

Y vT f 

11 5 ‘ 


Ygtanh(w-jahi )h2(a) + Y-iQtanhCwgah-j )h-j (a) 

L 

Y9tanh(w,ah, >1112(0) + v,otanh(W3ah, )m, (o) 


Iw, 


Ygtanh(w^ah-j )n2{a) 


I 


Y^Qtanh(w3ah-|)n^(a) 

Y-jQtanhCwgah^ )g-j (a) 


Y-j Qtanh(w3ah-j ) 


N 

^^.(h,-n)a/F?/lw3l 

Y , 

J r 

1 z{ 


Ygtanh(wiah-|)g2(a) 


- sign(w3)$^0 


12 


Ygtanh(w^ah-| )h2(a) 


h-,(a) 


- 


Y 

» 

Y„ 

1 1 

2 

^12 

U 3 I 


jtanh(w^ah^ )m2(a) + Y-| Qtanh(w3ah-j )m^ (a) 


ORIGINAL 
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+ ^Ygtanh(w^ah-j)n2(a) 


sina. 


TW, 


da 


" o * Jo 


^^13 



,-(h2-n)a/^/ |wtl_ -(h2+n)ot' 


>^/i«ti} 


• f 

- YQtanh{w,ah^)g2(ct)+Y^0^®^^^'^3“^l^^l^“^ 
. ^ 


sign(w-|)p 9 * 


Ygtanh(w-jah-j )h 2 (ot) + Y-io^anhCwgahi )h-| (a) 


* 

y ^ f 

'll' 


■ Iwtl 


Ygtanh(w-|ahi ) 


Y* ^ 6 ^ ^ 

n(i2(a) + ^|^| ~ YgtanhCw-jah^ )n2(ot) 




- e 


■ (h 2 +n)o‘v^/lw 3 || I ^ 


y 


-{h 2 -Ti)a/^/lw 3 


Ygtanh(w^ah^ ) 92 (ct) 


12 


+ Yigtanh{w2ahi )g-| (a) 

sign(w 3 ) 3 *Q Ygtanhlw^ah^ )h 2 (ct) + 

I ^ 

■' 

' 

Ygtanh(w^ah^ )ni 2 {a) + Yi 


III 

Y *2 


+ 


Yti y:x, 




Y X 

+ 11 Ps ■ ^ 


w. 


Ygtanh(w-jah^ )n 2 (a) + Y-j gtanhCw^ah^ )n-| (a) 


sina; 


da 
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APPENDIX C 


Evaluation of some Integrals 


e”^^sin bt dt = 


f a > 0 


••St A 

e cos bt dt = , a > 0 


x;;r cos ax dx = ^ 




roo 


0 rV »* <'* ' 7 




0 


rzuzr Sin 


ax dx = (1 - e"®^) 


23 


Iq e'^^cosh ay da = , a>0 
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APPENDIX D 


Derivation of the Stress Intensity Factors 
D.l Case of Integral Cracks ; 


The stress intensity factors are defined as: 


= Jim /2{x^-a) a^y(xpO) 


x.j->^ a 


kj^ “ liin / 2 ( X2“b) cf2y{x2»0) 


X2-^b 


(D. la ,b) 


From Eq. (4.16a) we can write: 
2 y^ >. ra 
V'yx' 


■ iiO-VxyJ ) 1, “ly(*r°> 


where o°y(xi,0) is a bounded function. 


Tri/2 


4,(t) = -Elii ruk 

(t-a)**{t+a)^ 


(D.2) 


Define th.j sectionally holomorphic function, 


i>{z) = 


IT 




iLlI 

t-z 


dt 


From [14] we obtain: 

Mz) = 


TTi/2 


, , - — i-+ iJ (z) 

(2a)^(z+a)^ (2a)^(z-a)^ ° 


(D.3) 


Using (D.3) to evaluate (D.2) and with the definitions 
given in (D.la,b) we have: 


-no- 



k - J- F(a) - 

* ■ ■ "VV' ’ 


and similarly 


k = . i F*(b) 

D (l_v* V* ) /b 
xy yx' 


' xy'^yx' 


★ * 

Fn(l) 

(1-V* V* ) 0 

' xy yx' 


D.2 Case of Broken Laminates 


'‘a = x’i^2 ^^‘*2 * >’2>^ “2y(*2-°> 

From Eq. (5.2b) 

2y* E* h 

a«^(x„,0) = Jl - j[ ' — f ^kpi.(xp,t)({)(t)dt + Op (xp.O) 

' (i-Vyx’ -'’1 

where o2 (xg*^) is a bounded function. 

J 


(D.4) 


(D.5) 


Define: 


<{)(t) = 


F(t) 


Ilili 


TTiy 




then : 


^(Z) = 1 


f ^ ^ dt 
Jk. t-z 

-hi 


ij^{z) = 


F(-h,)e”^ 1 


F(h,) 


{2hj )^s1nny (z+h^ (2h^ )^s1nirY (z-hj ° 


L— + 4 (z) (D.6) 


Using (0.4), (0.5) and (0.6) we obtain: 

* * 
lh,)'l ■ 

k 


(hi)V,4E/„(1) 


Wi1//87 


Wi 


1//^ 


'101 


Iwtl IwJ 

Y n02 ■ 

iwlllwi! 
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+ X 


103 


I“.|/ve7 , 

A. \ 

IM//S7' 

IWailwtl 

Y n04 

■|Wsllw:i 

Y 

J 


j 



kjj has the same expression as the one we derived in the 
previous case. 

D. 3 Case of a Crack Crossing the Interface 

D.3.a Singular Behavior at the Interface : 

In this case the definitions of the, stress intensity 
factors are: 


and 


kb = a2y(*2.0) 

S, 


‘xx = 


'‘xy= ^SxyC’rJ'’ 

From Eq. (6.2b) we can write: 

a2y(X2,0) = f ^ 

(1-v V )tt -'c > ^ 

^ xy yx^ 


(D. 7a ,b,c) 


(D.8) 


where a 2 y(x 2 » 0 ) is a bounded function. 




(z) 


(h 2 -c)®simr /2 (z-c)*® 


(h.-cj’^sInwB (z-h,)® ^ 


(D.9) 


Using (D.9) and (D.8), Eq. (0.7a) yields to: 

k. = F!i:ci^ 

(h 2 -c)^ 

or 

^ EjChj.c)*® . 

*^b ° V — Fj-i) 

<’-Vyx> 

Separating the singular parts of Eqs. (e.lSa.b) we have 

rh. 


X ih_ 


1 


j^^K^sCn.y)*^ (n)dn + <^°j^(y) 


2 G^y ^^xy^^^y) = 


rh 




K 3 s(n.y)<J)(Ti)dn 


f^2 * 

1 ^ K4s(n.y)4> (n)dn + T°^y(y) (D.lOa.b) 


where cr^^(y) and T°^y(y) are bounded. 


Let; 


■f(t) “ 


KlI 


F(t)e 


iTTg 
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then 



F(-h,.)e*^g , 

(2h^)®sinn6 (z+b^)® 


F(h,) 

(2h^)^sinTTe(z-h^)^ 


+ 'l'o(z) 


(D.ll) 


Using (D.ll), (D.9), (D.lOa.b) and the definitions (D.7b,c) 
after lengthy algebra we obtain: 



and 
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i 



"b = <’2y(X2.0) 

The derivation is similar to the one done for the case of 
internal cracks. Thus, we have: 



•2y*4E*^ 

( 1 -V V ) 

^ xyyx' 


V’> 
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